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We discuss a microscopic mechanism of the spin-triplet superconductivity in the quasi-two- 
dimensional ruthenium oxide Sr2Ru04 on the basis of two-dimensional three-band Hubbard 
model. We solve the linearized Eliashberg equation by taking into account the full momentum- 
frequency dependence of the order parameter for the spin-triplet and the spin-singlet states, 
and estimate the transition temperature as a function of the Coulomb integrals. The effective 
pairing interaction is expanded perturbatively with respect to the Coulomb interaction at the 
Ru sites up to the third order. As a result, we show that the spin-triplet p-wave state is 
more stable than the spin-singlet d-wave state for moderately strong Coulomb interaction. Our 
results suggest that one of the three bands, 7, plays a dominant role in the superconducting 
transition, and the pairing on the other two bands(Q! and /?) is induced passively through 
the inter-orbit couplings. The most significant momentum dependence for the p-wave pairing 
originates from the vertex correction terms, while the incommensurate antiferromagnetic spin 
fluctuations, which are observed in inelastic neutron scattering experiments, are expected to 
disturb the p-wave pairing by enhancing the d-wave pairing. Therefore we can regard the spin- 
triplet superconductivity in Sr2Ru04 as one of the natural results of the electron correlations, 
but cannot consider as a result of some strong magnetic fluctuations. We will also mention the 
normal Fermi liquid properties of Sr2Ru04. 

KEYWORDS: electron correlations, three-band Hubbard model, third order perturbation theory, vertex correc- 
tions, Sr2Ru04, spin-triplet superconductivity, superconducting transition temperature. 

§1. Introduction 

Clarifying the mechanism of the spin-triplet superconductivity in Sr2Ru04 has been one of the 
most challenging issues in the physics of strongly correlated electrons This unconventional 

pairing state has been accepted widely as a result of many intensive experimental and theoretical 
studies. 

Some of the most remarkable results supporting the odd-parity parallel-spin pairing were obtained 
by the direct measurements of the spin susceptibility through the superconducting transition. The 
NMR Knight shifti^ and the inelastic polarized neutronll' scattering measurements revealed that 
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the spin susceptibility of the electron system does not show any change at the transition. These 
facts indicate that the electron system does not lose all the activity in the spin space even in 
the superconducting phase, and exclude the possibility of the singlet Cooper pairing. In addition, 
there are many other observations suggesting indirectly but very strongly the realization of the 
unconventional pairing symmetry. The superconductivity is crucially destroyed by the slight con- 
centration of nonmagnetic impurities,& in contrast to the conventional s-wave superconductivity. 
The NMR relaxation rate exhibits no coherence peak just below The muon spin relaxation rate 
is increased below T^., which indicates that an internal magnetic field is turned on spontaneously 
and the time-reversal symmetry breaks down in the superconducting state.0^ More recently there 
are several experiments reporting that the second superconducting phase transition occurs in the 
vicinity of curve 

The normal state properties of Sr2Ru04 are quite different from those of high-Tc cuprates, in spite 
of the similarity in the crystal structures. Sr2Ru04 possesses a metallic conductivity without any 
carrier doping, and exhibits typical quasi-two-dimensional Fermi liquid behaviors, for example, the 
Pauli paramagnetism, the behavior in the resistivity at low temperaturesJll0' The electronic 
structure is characterized by the three cylindrical Fermi surfaces, a, (3 and 7, according to the 
de Haas-van Alphen measurements. 0) This is in good agreement with the results of the band 
structure calculations. 

Prom the theoretical point of view, Rice and Sigrist predicted the parallel-spin pairing shortly af- 
ter the discovery of the superconductivity0 Their discussions are based on the facts that Sr2Ru04 
gives a set of Fermi liquid parameters comparable to those of the superfluid ^He and that the three- 
dimensional version SrRuOs is a ferromagnet. According to that situation, it seems quite natural 
to expect that some ferromagnetic spin fluctuations induce the spin-triplet superconductivity. 
However the inelastic neutron scattering measurements did not detect any enhancement of ferro- 
magnetic components but observed a sizeable incommensurate antiferromagnetic fluctuations in 
the spin susceptibility x(9)-0^ Therefore it is difficult to expect that some strong ferromagnetic 
fluctuations play the leading role for the superconductivity of Sr2Ru04. 

Here we would like to comment on some other theoretical considerations on the mechanism of the 
superconductivity. Some theorists consider that the observed incommensurate antiferromagnetic 
fluctuations mediate the superconductivity. According to the discussions, the incommensu- 
rate fluctuations should possess a large anisotropy, XzziQmc) > 2% \-{Qmc)- To our knowledge, 

however, such a high anisotropy has never been observed in Sr2Ru04. Particularly, the recent 
neutron scattering measurement by Servant et al. suggests that the spin susceptibility x(Qmc) 
is very isotropic© Recently magnetic properties of Ti-doped Sr2Rui_a;Tij:04 have been investi- 
gated.!!!!) This compound exhibits the long range magnetic ordering for x > 0.03, which evolves 
from the incommensurate antiferromagnetic fluctuations. In the ordered phase, the spin moment 
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points along the c-axis directionES^ This anisotropy is quahtatively identical to that required in the 
anisotropic spin fluctuation scenario. What we should note here is, however, the following point. 
The effect of Ti-doping on the superconducting transition temperature is quite similar to that of 
other nonmagnetic impurity doping, in spite that the incommensurate fluctuations are much en- 
hanced by Ti-doping.il' If the incommensurate fluctuations mediated the superconductivity, the 
transition temperature should behave for Ti-dopings in a different way from the case of other non- 
magnetic impurity doping. Consequently, we consider that the incommensurate antiferromagntetic 
fluctuations are not essential for the superconductivity of Sr2Ru04. 

A scenario based on the incommensurate orbital fluctuations was discussed. In order to explain 
the spin-triplet superconductivity, however, they require rather strong inter-orbit repulsion between 
Ru4de electrons. The effect of the inter-orbit repulsion usually becomes weak remarkably by the 
screening. Therefore we consider that the incommensurate orbital fluctuations are not large enough 
to stabilize the spin-triplet pairing in Sr2Ru04. 

There are some other scenarios based on the Hund's coup ling .0' 00* The Hund's coupling 
is usually attractive only for two parallel-spin particles which are in the different orbits on the 
same localized atom. Therefore it seems that the Hund's coupling does not necessarily play an 
important part for stabilizing the parallel spins over the coherence length. In addition, the hy- 
bridization among the atomic orbitals Ru4de is relatively small in Sr2Ru04, because of both the 
two-dimensionality and the symmetry of the atomic wave functions. We consider that the Hund's 
coupling is actually hard to be the main origin of the spin-triplet superconductivity in Sr2Ru04. 

Recently the present authors discussed on the basis of the one-band Hubbard model for the single 
branch 7.® According to the previous works,!!® the spin-triplet superconductivity of Sr2Ru04 
is a natural result of the electron correlations. Even if the spin susceptibility has no prominent 
peak around q = 0, the effective pairing interaction has a momentum dependence favorable for 
the p-wave pairing state. The aim of the present article is to extend the previous single-band 
discussion to three-band one and give a comprehensive understanding on the mechanism of the p- 
wave pairing. As a result, we show the following points: the spin-triplet p-wave state is more stable 
than the spin-singlet d-wave state for moderately strong Coulomb interaction. One of the three 
bands, 7, plays the dominant role in the superconducting transition, and the pairing on the other 
two bands(a and /3) is induced passively through the inter-orbit couplings. The most significant 
momentum dependence for the p-wave pairing originates from the vertex correction terms, and 
the incommensurate antiferromagnetic spin fluctuations are expected to disturb the p-wave pairing 
by enhancing the d-wave pairing. Therefore we can regard the spin-triplet superconductivity in 
Sr2Ru04 as one of the natural results of electron correlations, and cannot consider as a result of 
some strong magnetic fluctuations. 

This paper is organized as follows. In § we give the formulation including introduction of the 
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Hubbard Hamiltonian and perturbation expansion. In § ^, we show the results of the numerical 
calculations of Tc as a function of the Coulomb energy, anomalous self-energy and effective inter- 
action. There we mention also normal Fermi liquid properties, Fermi surface, normal self-energy 
and density of states. Finally, in § ^, we comment on the results, give a proposal for experiments, 
and conclude the paper. 

§2. Formulation 

2.1 Three-band Hubbard model 

Sr2Ru04 is one of typical strongly correlated electron systems where Ru4d electrons play the 
most significant role for the electronic properties. According to the band calculation results, Ru4de 
orbitals take the main part of the density of states near the Fermi level, although the Ru4de and 
02p orbitals are hybridizing. EE) Therefore the electronic structure is expected to be reproduced 
well by considering three Wannier states, xy, yz and xz, which have the characters of Ru'idxy,yz,xz- 

We give the Hamiltonian in the form, 

H = Ho + H', (1) 
where Hq and H' are the noninter acting part and interaction part, respectively. Hq is given by 

where c^(c) is electron creation(annihilation) operator, and k, a^') and I denote momentum, spin 
and the Wannier states {xy,yz,xz), respectively. We take the following band dispersions, 

ixy{k) = 2ti{coskx + cosky) + At2CoskxCosky — fi^y, (3) 

^yzik) = 2t3COsky + 2t4COsA;a; - fiyz, (4) 

(.xzik) = 2t3Coskx + 2tiCosky - fixz, (5) 
A(fc) = 4Aosin/c2;sin/i:j,. (6) 

ti, t2, ts, t4 and Ao are transfer integrals between nearest or next nearest neighbor Ru sites (Fig. 1), 
where we have assumed tetragonal symmetry, and fi£ is the chemical potential. For the present 
calculation, we take h = -1.00, ^2 = -0.400, is = -1.25, ^4 = -0.125 and Aq = -0.200, and 
determine the chemical potentials fi£ to satisfy Uxy = 0.700, Uyz = 0.700 and rixz = 0.700, where 
ni is the electron filling per one spin state in the orbital £. The total population of the conduction 
electrons at a Ru site equals 2{nxy + Uy^ + nxz) = 4.200, which is a little larger than the electron 
number obtained by the experiments,© 4.032. We discuss this deviation in § ^. We take the 
following on-site Coulomb interaction for the interacting part. 
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i e^e' cr,cr' i i^fj cr^cr' 

J is the Hund's coupling. For the convenience of calculation, we express the Coulomb integrals by 
a matrix -^(^io-i)(£2o-2),(^30-3)(^4<t4)! where the matrix / has the following elements, 

h 

h 
In 



{ecr){ea),{ea){ecr) 


= u 


(fo)(£'a),{£'cr)(£cr) 


= U'-J 


{ecr){e'a),{e'a){ecr) 


= U' 


(fc)(£'CT),(£CT)(£'cr) 


= J 


(fc)(£a),(£'CT)(£'<T) 


= J' 



(8) 



where i ^ i', and a ^ a, and the other elements are zero. Using the matrix /, H' is expressed in 
the simple form, 

^'=2^ Y ^CiC2,C3C4 4ci4c2Cif3CiC4> (9) 

i C1C2C3C4 

where we have adopted the short notations, ( = {la) and '}2c^ = Z^^So-- We define the antisym- 
metric bare vertex by 

^ClC2,C3C4 = -^ClC2,C3a - -^ClC2,C4C3- (10) 

In the perturbation expansion, this vertex is represented diagrammatically by the empty square 
as depicted in Fig. 2. If we have in mind the rotation invar iance of the de atomic orbital wave 
functions in orbital space, the Coulomb integrals, {/, U' , J and J', are not independent parameters, 
and the relations U = U' + 2J and J = J' are derived.© We cannot, however, easily tell that 
the relation should be valid even for the Wannier wave functions in metallic solids, because the 
Coulomb interaction may be changed effectively, in general, for example, due to the screening effect. 
In the present calculation, therefore we treat the Coulomb integrals as independent parameters. 
In any case, as we can see in § |3.3| , we obtain the p-wave superconductivity for a broad parameter 
region of the inter-orbit couplings, and may conclude that this assumption of the independence is 
not significant in the present discussions. 

In the previous work, we investigated the magnetic properties of layered ruthenates on the basis 
of the same model.i^) 

2.2 Perturbation expansions 

We expand the normal self-energy and the effective pairing interaction perturbatively up to the 
third order with respect to the Coulomb interaction H' . At first, we neglect the hybridizing term 
including \{k). In Sr2Ru04, \{k) is the second nearest hopping term and expected to be small 
compared with the band width. By this simplification, the normal self-energy and the effective 
interaction are expanded perturbatively by the bare Green's functions, 

Gf\k) = -. (11) 
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where i = {xy,yz,xz} and k = (fc,ia;„). 

In the present scheme, the normal self-energy has only the diagonal elements The pertur- 

bation terms to be summed up are displayed in Fig. 3. Since the terms represented by Fig. 3(N1) 
only shift the chemical potentials /x^, we regard them as already included in /x^'s. The general forms 
of perturbation terms remained for summations are given as follows. 

2 



vN2 



\ J fel,fe2 C3C4 €l€2^3?4 

sF.^ctw = -(§)' E E E E Gg)3(^i)Ggk(^^)41(^-^i + ^^) 

\ / A:i,A:2,fc3C3C4?l?253?4 71727374 

xG(l(fc3)G(l(fc - /.I + fc3)r[S2,3C4rgU3e4rS3C.74' (i^) 



\ / *;i,fe2,fc3C3C4 66654 71727374 
xG(l(fe3)GS),,(fe + fei - fc3)rg),,3,/g),,3,/(°),,,3,,. (14) 

The factors \ and \ multiplied in the terms ^^^^^(fc) and are necessary to avoid the 

redundancy in the summation. By the simplification mentioned above, 

Gg)(A;) = Gf (15) 

where C, = {£cr), (' = {i'a'), and (5^^' = Sw^aa'- The diagonal elements of the normal self-energy 
are reduced to 

^ iV EE^y^^l) 2 E ^miU2(^ ~ E E^(M(m2i^2),(miVi)(£2<72)^(mii/i)(£2a2),(M("l2l'2 

E E -^muii^ ~ ^l)^mim2(^ ~ ^1) X] E E^(°a)(n2r2),(niri)(^20-2) 



nin2 mim2 T1T2 1/11^2 cr2 

(niTi)(m2!^2),(mii/i)(n2T2) {miUi){e2<T2),{icr){m2i^2) 

-i E E 4%2(^ + fci)^^U2(^ + ^1) E E ErSS)(,,.,),(„,.,){n2^) 

nin2 m\m2 T1T2 vii/2 <T2 

(niTi)(n2T2),(fnii^i)(fn2t'2) (mii^i)(m2i^2),(M{^2 0-2) 



(16) 



where the functions, xfj{q) and ^^^}{q), are calculated by 

4°' (^) = -^rT.Gf\Q + k)Gf\k), (17) 



*S'^(5) = -^EGr(^-^)G'?^W- (18) 



k 
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The renormalized Green's functions are given by 



G,{k) = -. (19) 

where = ^^(fe) + '^t{k) — Sfi^. S/ie is determined to satisfy the relation 

Sne = ^Y.(Geik)-Gf\k))=0. (20) 

k 

Here we take into account the effect of the hybridizing term A(fc). We consider the following 
renormalized Green's functions: 

Gxy,xy{k) = T— 1=^ 77T, (21) 

^y^'V'^^^ = : ^-7: A^(fc) ' ^^^^ 

Gyz,xz{k) = Gxz,yz{k) (24) 

^ m ^ 

(iWn - ^yz{k)){^n - ^xz{k)) - A2(fc) ' 

and the other elements of Gui{k) are zero. These Green's functions are the matrix elements of the 
operator [iLVn — Hq — {Ti — 5^)]~^ , where S — is an operator with the diagonal matrix elements 
(E£(fe) — 6ni)5w. Diagonalizing this renormalized Green's function, we obtain the following form 

G,{k) = : (25) 

Here v = {a,/?, 7}, and 

E^,f,{k) = E+{k) T mk) + A2(fe)]i (26) 

E±ik)^^iEyAk)±Exzik)) (27) 

E^{k) = Exy{k) (28) 

Then we expand the anomalous self-energy perturbatively in a similar manner. The perturbation 
terms to be summed up are displayed diagrammatically in Fig. 4. The general forms for the 
perturbation terms are given as follows. 

4cIw = -^^EE4c3(^')r2L,c.a' (29) 

k' C3C4 

^acI(^) = -f§)'EE E 4cs(k')G%{k^)G(Xik + k^-k')^^^^^^^ (30) 

\ / fe'.fci C3C4 71727374 
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\ / A;',A:i,fe2C3C4a6?354 71727374 



^^7472 (^2)G^3li {k k' + A;2)rl I J , , I ^ (31) 



X<-^7472 V"':^;^737l V"- '^■2;^ C46,S3Cl'- €l72,73«4'- 7lC3,C274' 



^c^c?(^) = f§)' E E E E 4c3mi(-^ + ^' + ^i)GS^74(^i) 

\ / fe',fcl,fe2C3C4ei6^3?4 71727374 

x^SLl-^ + + fc2)G(2^(fe2)r2L73Ci<k,C2«342C3,C47i' (32) 



^fic?(^) = f^)' E E E E 4c3mi(fc-fe' + fci)GS).(fci) 

V / fc,fel,fe2C3C4 66C3C4 71727374 

xGg)^(fc + k, + fc2)Gg)^(fc2)r(l,C3CirS3,c27irgL,a73' (33) 
^Cicf(^) = K^)' E E E E 4c3(^')Ga(^-fc' + ^i)Gil(^-i) 

\ / A;,fei,fc2C3C4 666^4 71727374 

xGg3(fc + k,- ^2)41(fc2)r[l,,,,^rS3,,,,2rgL,..c2' (34) 



^ci6 (^) = K§) ' E E E E 4c3 (^OG^) 3 (-fc + ^' + k.)G^X (ki) 

\ / fc,fcl,fe2C3C4€l6€3€4 71727374 

xG^li-k + k,- ^2)^^)3 (fc2)rg),,,,,,r[f,3,,,,,<^,^^^^, (35) 

where F^^^(/c) is the anomalous Green's function. The factors ^ multiphed in the (Al), (A3d) and 
(A3e) are necessary to avoid the redundancy in the summation. By summing up these contributions, 
we obtain the following Eliashberg equation within the third order perturbation theory. 

^fil2(^) = -^E E 4c3(^')(rc4C3.C2Ci(^',^) - ^rSkc2Ci)' (36) 

k' C3C4 

where T^^^g i^jCi (^'' k) is the renormalized pair scattering amplitude and given in the following, 

rC4C3,C2Cl(^'>fc) = r[4[3,^2Cl ~ E ^mlmai^ " ^')^C473,7lCl^7lC3,C273 

7173 

+ E E ^niUii^ - k')X!^\m2ik - k')^[Ai2,e.lC,l^\il2,lli2^l}(:3,C,212 
7172 ^1^2 



lT.m:Gt\{-k+k'+k,) 

ki 7173 66 



^n2ni { «; «1 C473,7iCi 7i6,C26^ 6C3,673 
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+ ^'^n2nA 1^ + f^l)^ 66,Ci7i f473,65i 7iC3,C273 



T 



1 



ki 7173^1^2 
,{0) 



.(0) 



.(0) 



73C3,f271 C2C4,5l73 



(0) 



p(0) 

2-™2'"lV- ' •--^'^ C47l:73Cl 73C3,66 Cl6>7lC2 



.(0) 



(37) 



Here we have adopted the short notations, Q = (iiCi), ji = {rriiUi) and = (riiTi). We assume the 
pairing interaction binds the most strongly the quasi-particles on the same band, and consider the 
three components of the anomalous Green's function and the anomalous self-energy, F^^^, (k) and 
T.^}^,{k), where = q,/3,7. The diagonalized band indices z^(=a,/?,7) are replaced by the orbital 
indices £{=xy,yz,xz) through the following relations. 



= E u\yk)ufl{k)\G^{k)\^^±,{k) 



(38) 



(39) 



where C = (^cr), C = {t'^^'), and JJ^^\k) is a matrix for the diagonalization of the noninteracting 
Hamiltonian i^o, 

7 a (3 



where 



yz 
xz 



1 
K{k) L{k) 
-L(fc) K{k) 



(40) 



M(fc) 



L(fe) = sgn(A(fe))\/-(l + M(fe)), 



e(fc) + A2(fc)- 



.At 



(41) 
(42) 
(43) 

(44) 



We replace the left hand side of the eq. ( p6[ ) by A(T)S^^^2 (g), and determine the transition temper- 
ature, at which the eigenvalue A(T) equals unity (A(Tc) = 1.00), by solving the set of the equations 
eqs. ([3^)-(39). It is convenient to express the anomalous self-energy in the following form© For 
triplet pair, 

i:t.'ik) = KD,{k).a)<jyU,, (45) 
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and for singlet pair, 

j:ta'ik) = [i^Ak)cTyU', (46) 

where cr is the Pauli matrix, and the function, Du{k) {'^i,{k))^ behaves as a vector (scalar) un- 
der the rotation in the spin space. In the numerical calculation, we assume the vector Du{k) is 
perpendicular to the basal plane, 

D^{k) = D^{k)z. (47) 

The direction is consistent with the experimental results of Knight shift.!) This assumption is 
not substantial in the present calculation, because we do not introduce any anisotropy in the spin 
space. We can show that, if we take other direction of the vector, we obtain the same transition 
temperature. In order to determine the direction of the vector Diy{k), we must take into account 
the effect of the spin-orbit coupling. We do not consider that the spin-orbit coupling is essential 
for discussing a mechanism of the spin-triplet superconductivity in Sr2Ru04. 
Here we define the function, T^'^_^(^^^,y^^^-^^{k' ,k), by 

As we see in § the component for the parallel spin pair on the 7 band, r^o-(T,7(T(T(^'', k), possesses 
a momentum dependence favoring p-wave pairing. 

In conclusion of this section, note that the normal self-energy and the effective interaction have 
been expanded up to the third order in Coulomb energy by using only the bare G^^\k) {£ = 
{xy, yz, xz}), but the quasi-particles forming Cooper pair are described by the renormalized Gy{k) 
V = {a,/5,7}. 

§3. Results of Calculation 

3. 1 Details of numerical calculations 

The momentum and frequency summations appearing in the equations, (16)-(18), (^) and (37) 
can be performed numerically with use of fast Fourier transformation algorithm. For the numer- 
ical calculations, the first Brillouin zone is divided into 128 x 128 fc-meshes, and the number of 
Matsubara frequencies taken is A'^^ = 1024. The temperature region for numerical calculations is 
bounded from below by T > W/2'KNf ~ 0.0012, where W is the noninteracting band width. All of 
the numerical calculations in the present paper are performed in T > 0.00250. This lower bound 
corresponds to about 10(K), if we assume the band width W is about 2(eV), and is still higher 
than the real of Sr2Ru04, 1.5(K). Therefore we must extrapolate the calculation results of 
to the weak interaction region. 

3.2 Normal Fermi liquid state properties 

We show here the results of the quantities related to the normal Fermi liquid state. 
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The Fermi surface consists of three sheets as depicted in Fig. 5. The circle around the point 
(7r,7r), a, is hole-like, while the large circles around the corners, /? and 7, enclose electrons. These 
results are in good agreement with those obtained by de Haas-van Alphen measurements0^ and 
recent angle-resolved-photoemission measurements. 

The normal self-energy on the 7 band as a function of frequency, 'E,^{kp,ij), is shown in Fig. 6. 
We can see usual Fermi liquid behaviors for every case of various inter-orbit couplings, ReS(LL') ~ oj 
and ImS(u;) ~ c'w^ (c, c' < 0) at low energy (uj ~ 0). The remarkable property which we would 
like to point out here is that inter-orbit couplings do not affect the normal self-energy in the low 
energy region u; ~ 0, while they do in the high energy region. This is partly due to the situation 
where the hybridization among the localized orbits, xy, yz and xz, is not so large in Sr2Ru04. As 
far as we discuss low energy phenomena, such as superconductivity, therefore, we may expect that 
the single-band approach@ is a good starting point. 

Next we show the density of states(DOS) for a typical case in Fig. 7. The DOS of each band u 
(= a, (3, 7) is calculated by the formula, 

p,{u;) = --Y,lmG^{k,u;), (49) 

where G^{k) is obtained by continuing analytically Gu{k) from the upper half plane. What we 
should note is the large DOS peak near the Fermi level. This peak originates from the van Hove 
singularity on the two-dimensional band 7. As pointed out by Singhjll' if the Fermi level comes 
above the singularity, the Fermi surface 7 changes its topology, from electron-like circle to hole-like 
circle around the point (7r,7r). The 7 band takes the main part of the DOS at the Fermi level. This 
agrees with the results of the de Haas-van Alphen measurements. As a result, the 7 band is 
expected to be significant in discussing the electronic properties of Sr2Ru04. 

3.3 Eigenvalues for triplet p-wave and singlet d-wave states 

We show the calculation results of the maximum eigenvalues for the spin-triplet and the singlet 
states in Fig. 8. There we show the calculation results for four cases of relatively strong inter-orbit 
couplings. The orbital symmetry of the state giving the maximum eigenvalues is p-wave for the 
spin-triplet state and d^2_y2-wave for the spin-singlet state. In the cases (a), (b) and (d), we obtain 
the larger eigenvalues for the spin-triplet state than for the spin-singlet state at low temperatures. 
Therefore we expect that the p-wave state is realized for the three cases. For the case (c), where the 
Hund's coupling is strong, we can see that the eigenvalues for the p-wave state is small compared 
with those for the d^2_y2 state. The reason why the d^2_y2-wave pairing state becomes more stable 
than p-wave one, as the strength of the inter-orbit coupling increases, will be discussed in § It 
should be noted here that the values of inter-orbit couplings in the cases, (a), (b) and (c) are too 
large. We consider that actually they should be smaller, since in general they are weakened more 
by the screening effects than the intra-orbit coupling U. If we take smaller realistic inter-orbit 
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couplings, for example, as in the case (d), (i^2_j^2-wave state is not expected to be realized at all. 
Note that the situations for which we have given the results of eigenvalues here are still relatively 
unfavorable for the p-wave state because of the strong inter-orbit couplings. 

3.4 Transition temperature for triplet p-wave state 

The calculated transition temperature for the p-wave state is displayed in Fig. 9 as a function 
of U. In the both cases of the inter-orbit couplings, (a) and (b), we may obtain higher transition 
temperature for the p-wave state, according to the results of calculation of eigenvalues in § |3.3| . 
We consider that the lower bound of temperature for the calculation, T = 0.00300, corresponds to 
about 10(K). We may obtain the more realistic value of Tc=1.5(K), by extrapolating the results 
down to the smaller U region. According to the results, Tc becomes high for strong inter-orbit 
couplings. In the limit of weak inter-orbit coupling, we obtain the same results as we did in the 
previous single-band analysis for the 7 band. 

3.5 Momentum dependence of anomalous self-energy parts 

The momentum dependence of the anomalous self-energy Dy{k) is plotted in the Fig. 10, where 
we consider the state with the ky-\\ke symmetry Dy{kx,ky,\ujn) = —D^{kx,—ky,iu)n)- The results 
show highly anisotropic p-wave state, and not /-wave. If we assume in the beginning the condition, 
Di,(kx,ky,iuJn) = —Di^{—kx,ky,iuJn), we obtain the p-wave state with the kx-like symmetry. The 
relative phases of D^{k)^s converge to zero, and we can assume without any loss of the generality 
that D^{k)^s are real. 

We should note that the magnitude of the anomalous self-energy is the largest on the 7 Fermi 
surface. This is "orbital dependent superconductivity" , which was proposed originally by Agterberg 
et a/.il' The magnitudes of Da{k) and Dp[k) are almost the same, since both of the a and the (3 
bands are constructed by hybridizing the original one-dimensional yz and xz bands. The 7 band 
originates mainly from the two-dimensional xy band, and is separated from the other two bands, 
a and /?, in that the hybridization is expected to be small (In the present theoretical formulation, 
the hybridization is assumed to be exactly zero, since we consider only the single Ru02 layer). In 
Fig. 10, we have shown the case of relatively strong inter-orbit couplings. If we take small inter- 
orbit couplings, the relative magnitude of D^{k) to Dafi{k) becomes large. In the unrealistic and 
extreme case where inter-orbit couplings are zero {U' = J = J' = 0), Da{k) and Dp[k) converge 
to zero, and only D^{k) converges to a finite value. As a consequence, we find that the 7 band 
dominates the superconducting transition. Here we have shown only the results for the triplet 
p-wave state. According to our calculation, the dominance of the 7 band is also valid for the case 
of the singlet d^2_„2-wave state. 
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3. 6 Momentum dependence of effective pairing interaction on the 7 hand 

We extract only the pair scattering amplitude of quasi-particles from the 7 band to the 7 band. 
This channel is the most significant not only for the p-wave pairing but also for the d^i-yi-wave 
pairing. In Fig. 11 we show the plots of the vertex functions (a)r^o-(T,7(T(T(^') k) for the parallel spin 
pairs and (b)r^o-CT,7CTo-(^') ^) for the anti-parallel spin pairs. 

In the triplet case (Fig. 11(a)), we find that the function T^(j^^^(jfj{k' ,k), has a maximum value 
around the point k ~ —k'. This feature originates from the vertex correction terms (mainly, (A3d) 
and (A3e) in Fig. 4), and not from the exchange of any bosonic excitations. This characteristic 
momentum dependence of the vertex function is the origin of the relatively high transition temper- 
ature for the triplet pairing state, and is basically the same one as derived in the previous work for 
the single-band Hubbard model.il' 

In the singlet case (Fig. 11(b)), r^o-CT,7CTcr('?; k) shows the characteristic peak around k = A;' + Qiaf> 
where Qiaf ~ (±0.67r, ibO.Gvr). These peaks are due to the nesting of the quasi-one-dimensional 
Fermi surfaces, yz and xz (or almost equivalently, a and /?). This momentum dependence was 
observed as sizeable incommensurate antiferromagnetic spin fluctuations by the inelastic neutron 
scattering measurements, 

0) 

and was obtained also by the theoretical calculations on the magnetic 
properties of Sr2RuO4.00^ Since it is through the inter-orbit couplings that the momentum 
dependence of the incommensurate fluctuations is reflected in the pairing interaction on the 7 band, 
the peaks around k = k' + Qiaf vanish as the strength of the inter-orbit couplings is decreased. On 
the other hand, the large values of the vertex function around fc ~ fc' are due to the ferromagnetic 
components of the susceptibility for the xy{'y) band. Note that the values of the vertex function 
T-yf^^^^^f^{k' , k) are positive all over the Brillouin zone. This is because of the on-site repulsion U. In 
general this is the reason why conventional s-wave superconductivity is not expected to be realized 
in strongly correlated electron systems. 



Note that, as mentioned in § 3.3, the cases we have shown here are the strong inter-orbit coupling 
ones. We consider that the inter-orbit couplings are actually not so strong. We have assumed such 
strong inter-orbit couplings to elucidate their roles. 

§4. Discussion and Conclusion 

In this section, we would like to discuss the results of calculation comprehensively, and give some 
proposals on experiments, based on our microscopic picture of the superconductivity in Sr2Ru04. 

In § we have taken the electron filling per one spin state, n^y = Uyz = Uxz = 0.700. This 
seems to be a little large compared with the results of the de Haas-van Alphen measurements,0^ 
although the Fermi surfaces for this filling reproduce the observed ones qualitatively well, as shown 
in § 3.2. According to the measurements, it is plausible that the orbits are filled as the total 



electron number equals 4.032. If we assume that all of the three orbits are filled with even number 
of electrons, the electron number of one orbit is = 0.672. If we take the filling ui = 0.672 in the 
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calculation, the parameter region of the inter-orbit couplings where the p-wave state overcomes the 
d^2„y2-wave state is limited only in the weak inter-orbit couplings. We consider, however, that this 
situation is unrealistic, because (i^2_j^2-wave state has never been observed experimentally near the 
triplet superconducting phase of Sr2Ru04. We relax the relation 2{nxy + riyz + nxz) = 4.032 in 
order to obtain the good agreement with the realistic situation. The incommensurate nesting vector 
QiAF {={QiAFx^QiAFy, 0)) is estimated by Qiaf^ ~ 2(1 - nyz)^ and QiaFx ~ 2(1 - Uxz)-!^- The 
observed nesting vector is QiaFx = QiAFy = O.GOtt.EI^ Therefore we should take Uyz = rixz = 0.700, 
and W6 licivG assuniGci that the three orbits ctrc almost d-Cgencratc cLnd fixy — 0.700. At the present 
stage, the theory proposed here suggests that Sr2Ru04 seems to be located near the boundary of 
the p-wave and the dx2_y2 states in the parameter space. We consider, however, that the present 
theory may be still insufficient for discussing quantitatively the competition of the p-wave and the 
dx2_y2-wave pairing states. In the future, it must be proved by more refined theoretical treatments 
that the dx2_y2-wave pairing state is more suppressed. One way of the further study is to inspect 
the effect of the higer order terms. In any case, we believe that, as far as we discuss only the 
mechanism of the p-wave superconductivity in Sr2Ru04, the theory gives a satisfactory result, in 
that we have obtained the relatively high Tc for the p-wave state and the momentum dependence 
which favors the pairing state, as shown in Figs. 9 and 11(a). 

In § we have shown that for the strong inter-orbit couplings, particularly for the strong Hund's 
coupling, the eigenvalues for the dj.2_y2-wave state become large compared with those for the p- 
wave state. This is because the stronger the inter-orbit coupling becomes, the more prominently 
the incommensurate spin fluctuations are reflected in the pairing interaction on the 7 band. As 
the present authors pointed out in the previous work on the magnetic properties of quasi-two- 
dimensional ruthenates,iil' the Hund's coupling enhances the incommensurate fluctuations rather 
than the ferromagnetic components of the spin susceptibility. Therefore we may conclude that the 
Hund's coupling disturbs the p-wave pairing by enhancing the d-wave pairing. 



As shown in § p.5| , the orbital symmetry of the pairing is anisotropic p-wave. In general, if we 
assume the most promising form of the vector d{k) ~ {kxHky)z, 

i) 

we usually expect the nodeless 

gap around the cylindrical Fermi surface. The nodeless energy gap in Sr2Ru04 has, however, been 
considered to be inconsistent with the observed power-law behaviors in various quantities at low 
temperatures.i'i'ii'lil^ Recently we have successfully shown that the p-wave superconducting gap 
derived in the present formulation is quite consistent with the power-law behavior of the speciflc 
heat below Tc.© There the strong momentum dependence of the anomalous self-energy is essential 
for a node-like structure on the P Fermi surface, and results in the power-law behavior at the low 
temperature. Accordingly, we would like to point out here that the p-wave state obtained in our 
discussions can explain the power-law behaviors of the specific heat, even if we assume the symmetry 
d{k) {kx ± iky)z. 
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We would like to give some proposals for experiments, based on our microscopic picture on the 
superconductivity of Sr2Ru04. Let us consider that the number of electrons in the band ^{xy) is 
increased. In such a case, as we previously showed in the single-band analysis for the 7 band,® the 
superconducting transition temperature is expected to be enhanced. This is due to the situation 
where, if we increase the electrons in 7 band, the Fermi level becomes close to the van Hove 
singularity on the 7 band and the density of states at the Fermi level is increased. Note here that 
the van Hove singularity is located only slightly above the Fermi level, as shown in Fig. 7. At 
the same time, as we insisted in the previous work,© the ferromagnetic components of the spin 
susceptibility may be enhanced. Carrier dopings by chemical substitution usually not only change 
the carrier number but also damage the conducting RUO2 plane. Since the triplet superconductivity 
is sensitively destroyed by disorders, we should not take any chemical substitution for the carrier 
doping. We think that the most practical way of increasing the electron number in the 7 band 
without damaging the RUO2 layers at all is application of the uniaxial pressure along the c-axis. 
The orbits yz and xz spatially extends along the c-axis, while the orbit xy will be flat along the 
basal plane. If we apply the uniaxial pressure, the energy levels of the orbits yz and xz are lifted due 
to the crystalline field effect, while the xy orbit is not. Therefore we may expect that the electrons 
are transferred from the yz and the xz orbits to the xy orbit. Consequently, we would like to 
point out the following possibility. If we can indeed transfer the electrons from the yz and the xz 
orbits to the xy orbit by applying the uniaxial pressure, the transition temperature may become 
higher. At the same time, the system may become close to the ferromagnetism, as we suggested 
in the previous work.© Measuring the transition temperature and the uniform spin susceptibility 
under the uniaxial pressure along the c-axis is, therefore, very interesting. Recent experimental 
work on the elastic moduli suggests the possibility that the transition temperature is increased by 
applying the uniaxial pressure along the c-axis.ii^ We may consider the fact as an evidence that the 
superconducting transition occurs predominantly on the 7 band. This experimental fact supports 
the p-wave pairing rather than the d-wave pairing, because, if the (ia;2_y2-wave state was realized, 
the transition temperature should become lower by increasing the electron number of the 7 band, 
according to our previous work. 

In conclusion, we summarize the main results obtained. We have discussed a mechanism of the 
spin-triplet superconductivity in Sr2Ru04. We can obtain the p-wave pairing state for moderately 
strong inter-orbit couplings. There one of the three bands, 7, plays the dominant role in the 
superconducting transition, and the pairing on the other two bands(a and /3) is induced passively 
through the inter-orbit couplings. The most significant momentum dependence for the p-wave 
pairing originates from the vertex correction terms, and is basically the same one as we obtained in 
the previous single-band analysis© Therefore we should regard the spin-triplet superconductivity 
in Sr2Ru04 as one of the natural results of electron correlations, and cannot consider as a result 
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of some strong magnetic fluctuations. We would like to insist that Sr2Ru04 realizes the p-wave 
superconductivity which is basically described by a simple one-band repulsive Hubbard model. 
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Fig. 1. Schematic illustration of single Ru02 layer. The Ku4dxy,yz,xz and 02px,y,z orbitals construct the two- 
dimensional network by hybridizing with each other. 



Fig. 3. Perturbation terms for normal self-energy up to the third order. The shaded circle in (Nl) represents 
arbitrary self-energy insertion. The solid line and the empty square denote the bare Green's function G^^\k) and 
the bare antisymmetric vertex F^'^' , respectively. 

Fig. 4. Perturbation terms for anomalous self-energy up to the third order. The thin solid line, the thick solid line 
and the empty square denote the bare Green's function Gf\k), the anomalous Green's function F^^,{k) and the 
antisymmetric bare vertex F'"' , respectively. 

Fig. 5. The Fermi surface. The Coulomb integrals are U = 3.731, U' = 0.300(7, J = J' = U' , and the temperature 
is T = 0.0100. 

Fig. 6. The normal self-energy on the 7 Fermi surface as a function of frequency, (a) real part, (b) imaginary part. 

Fig. 7. The density of states. The Fermi level corresponds to oj = 0. (a) The total density of states. (b)The partial 
density of states for the three bands, a, /? and 7. The inset shows the details near the Fermi level. 

Fig. 8. The maximum eigenvalues as a function of temperature T. For all of the cases, the intra-orbit coupling 
is fixed ssU = 3.079. (a) U' = 0.500C/, J = J' = 0.667?/'. (b) U' = 0.500?7, J = J' = U' . (c) U' = 0.500Z7, 
J = J' = 1.33317'. (d) U' = 0.300Z7, J = J' = U'. In the cases, (a), (b) and (d), we expect that the spin-triplet 
p-wave state is stable compared with the spin-singlet dx2_y2 state. 

Fig. 9. Transition temperature as a function of U for spin-triplet p-wave state, (a) J = J' = 0.667J7'. (b) 



Fig. 10. Contourplots of the results of D^{k,iTTT) {i'=a, /?, 7). In the light and the dark colored regions, the 
functions take higher and lower values, respectively. The thick line circles represent the Fermi surfaces. The 
parameters are U = 3.385, U' = 0.500U, J = J' = U' and T = 0.00300. 

Fig. 11. Contourplots of effective interaction on the 7 band. The thick line circle around the corner is the Fermi 
surface 7. (a) For the quasi-particles with parallel spins, T^cra,-yrTrT{k' ,mT;k,mT). (b) For the quasi-particles with 
anti-parallel spins, Tjaa,-ysa{k' ,iT:T;k,mT). In both cases, fe' is fixed as pointed by the arrow. The parameters 
are U = 3.385, U' = 0.500[7, J = J' = U' and T = 0.00700. 



Fig. 2. The antisymmetric bare vertex. 
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Fig. 3. T. Nomura & K. Yamada 
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Fig. 4. T. Nomura & K. Yamada 




Fig. 5. T. Nomura & K. Yamada 
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Fig. 7. T. Nomura & K. Yamada 
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Fig. 8. T. Nomura & K. Yamada 
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Fig. 9. T. Nomura & K. Yamada 
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Fig. 11. T. Nomura & K. Yamada 



29 



